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Abstract. Motivated by the ubiquity of control-affine systems in optimal control theory, 
we investigate the geometry of point-affine control systems with metric structures in dimen- 
sions two and three. We compute local isometric invariants for point-affine distributions of 
constant type with metric structures for systems with 2 states and 1 control and systems 
with 3 states and 1 control, and use Pontryagin's maximum principle to find geodesic tra- 
jectories for homogeneous examples. Even in these low dimensions, the behavior of these 
systems is surprisingly rich and varied. 
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1 Introduction 

In [1], we investigated the local structure of point-affine distributions. A rank-s point-affine 
distribution on an re-dimensional manifold M is a sub-bundle 3" of the tangent bundle TM such 
that, for each x £ M, the fiber J x = T X M n 3" is an s-dimensional affine subspace of T X M 
that contains a distinguished point. In local coordinates, the points of 3" are parametrized 
by s + 1 pointwise independent smooth vector fields vq(x),vi(x), . . . ,v s (x) for which 3~ x = 
vq(x) + span (vi(x), . . . , v s (x)) and vq(x) is the distinguished point in 1 X . 

Our interest in point-affine distributions is motivated by a family of ordinary differential 
equations that occurs in control theory: the control-affine systems. A control system is a system 
of underdetermined ODEs 



where x £ M and u takes values in an s-dimensional manifold U. The system is control-affine if 
the right-hand side is affine linear in the control variables u, i.e., if the system locally has the 



x = f(x,u) 



form 



s 




(1.1) 



1=1 



where the controls u , . . . ,u s appear linearly in the right hand side and vq, . . . ,v s are s + 1 
independent vector fields (see, e.g., [3]). Replacing vq, which is called the drift vector field, with 
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a linear combination of vx,...,v s added to vq would yield an equivalent system of differential 
equations. In many instances, however, there is a distinguished null value for the controls (for 
example, consider turning off all motors on a boat drifting downstream), and this null value 
determines a distinguished drift vector field. In these instances, we always choose vq to be the 
distinguished drift vector field. Consequently, the null value for the controls will be 

u 1 = • • • = u s = 0. 

While the control-affine systems (1.1) may appear to be rather special, these systems are 
ubiquitous. In fact, any control system whatsoever becomes control-affine after a single pro- 
longation, so these systems actually encompass all control systems, at the cost of increasing the 
number of state variables. 

In [1] we studied local diffeomorphism invariants for these point-affine structures. A local 
equivalence for two point-affine structures is a local diffeomorphism of M whose derivative maps 
one distinguished drift vector field to the other, and maps one affine sub-bundle to the other 
(see [1] for precise definitions). With this notion of local equivalence, we were able to determine 
local normal forms for strictly affine, rank-1 point-affine structures of constant type when the 
manifold M had dimension 2 or 3. In some cases the normal forms are parametrized by arbitrary 
functions. 

The current paper seeks to refine the previous results by adding a metric structure to the 
point-affine structure. We do so by introducing a positive definite quadratic cost functional 
Q : 3" — > R. In local coordinates, where 

s 

w = v (x) + i>i(x)u J G %, 

we will define 

Qx{w) = ^2gij(x)u l u J , 

where the matrix (gij(x)) is positive definite and the components are smooth functions of x. 
This is a natural extension of the well-studied notion of a sub-Riemannian metric on a linear 
distribution, which represents a quadratic cost functional for a driftless system (see, e.g., [4, 5, 6]). 

With the added metric structure, we refine our notion of local point-affine equivalence to that 
of a local point-affine isometry. A local point-affine isometry is a local point-affine equivalence 
that additionally preserves the quadratic cost functional. 

Let 7(t) = x{t) be a trajectory for (1.1). The added metric structure allows us to assign the 
following energy cost functional to 7(i): 

E{l) = \J Q x{t) {x{t))dt. (1.2) 

Naturally associated to (1.2) is the optimal control problem of finding trajectories of (1.1) that 
minimize (1.2). We will use Pontryagin's maximum principle to find an ODE system on T*M 
with the property that any minimal cost trajectory for (1.1) must be the projection of some 
solution for the ODE system on T*M. 

In this paper we shall only consider homogeneous examples, i.e., examples that admit a sym- 
metry group which acts transitively on M. We shall use the normal forms from [1] as starting 
points, adding a homogeneous metric structure to the point-affine structure in each case. Even in 
these low-dimensional cases, the analysis can be quite involved; we will see that these structures 
exhibit surprisingly rich and varied behavior. 
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2 Normal forms for homogeneous cases 

We begin by identifying the homogeneous examples of the point-affine systems described in [1] 
where possible, and then we describe the homogeneous metric structures on these systems. In 
some cases, the metric structure must be added before the homogeneous examples can be iden- 
tified. Recall that the assumption of homogeneity is equivalent to the condition that all structure 
functions T l - k appearing in the structure equations for a canonical coframing are constants (see [2] 
for details). 



2.1 Two states and one control 

In [1], we found two local normal forms under point-affine equivalence. 

+ span (ggj 



Case 1.1. 3~ = + span (A), The framing 



d 1_8 

1 dx 1 ' V2 A dx 2 

(well-defined up to scaling in v 2 ) has dual coframing 

i] 1 = dx 1 , rj 2 = Xdx 2 , (2.1) 
with structure equations 

dn 1 = 0, drj 2 = mod rj 2 . 

Because the method of equivalence does not lead to a completely determined canonical cofra- 
ming, it is not clear from these structure equations whether this example is homogeneous as 
a point-affine distribution. 

Fortunately, this ambiguity is resolved when we add a metric function to the point-affine 
structure. This amounts to a choice of function G(x) > for which the quadratic cost functional 
is given by 

a a\ i_. , 2 



For the point-affine structure, the frame vector v 2 is only well-defined up to a scale factor; 
however, when we impose a metric structure (2.2), we can choose v 2 canonically (up to sign) by 
requiring that it be a unit vector for the metric. This choice leads to a canonical framing 

d 1 d 



dx^ z ^G{x)dx^ 
with corresponding canonical coframing 



r] 1 =dx 1 , if = ^G{x) dx 2 . 
The structure equations for this refined coframing are 

dri 1 =0, dr] 2 = I^t/ At? 2 , 

and so the structure is homogeneous if and only if -^r is equal to a constant c\. This condition 
implies that 

G{x\x 2 ) =G {x 2 )e 2c ' xl 

for some function Gq (x 2 ) . 
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The local coordinates in the coframing (2.1) are only determined up to transformations of 
the form 



x 



1 -x 1 + a, x 2 = cp(x 2 ), (2.3) 



and under this transformation we have 

G o (x 2 )=e 2 ^((cA'(x 2 )) 2 Go(0(x 2 )). 

Therefore, we can apply a transformation of the form (2.3) to arrange that Go(x 2 ) = 1, and 
hence G = e 2cix . Moreover, coordinates for which G has this form are uniquely determined up 
to a transformation of the form 



x 



1 = x 1 + a, x 2 = e~ cia x 2 + b. 



To summarize: the homogeneous metrics in this case are given by quadratic functionals of 
the form 

Q (JL+ U JL\ =I e 2 C1 *v 

\dx l dx 2 J 2 
for some constant c\, with corresponding canonical coframings 
r, 1 = dx 1 , i] 2 = e cixl dx 2 . 

Case 1.2. 3" = x 2 (^j^i + ^gfsj + span (gfs)- We found a canonical framing 

^ + J 7TT ' V 2 = x 2 - T , (2.4) 
ox 1 ox z I ox z 

with dual coframing 

r) 1 = ^dx\ rj 2 = ^(dx 2 - Jdx 1 ), (2.5) 
and structure equations 



dr] 1 =j) 1 Ai) 2 , di] 2 = T 2 2 r] 1 A r] 2 , 

where 

T 2 2 = x 2 lf 2 -J. (2.6) 

The structure is homogeneous if and only if T\ 2 is equal to a constant —jo- According to 
equation (2.6), this is the case if and only if 

J = x 2 J 1 (x 1 )+j (2.7) 

for some function J\ (x 1 ) . 

The local coordinates in the coframing (2.5) are only determined up to transformations of 
the form 



x = ft) I x 



), x 2 = xV^ 1 ), (2i 
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and under this transformation we have 

j(x\ x 2 ) = J(0(£ 1 ) , xV " * 2 ^l|- 

In the homogeneous case (2.7), this implies that 

4>"{x l ) 



J 1 {x 1 )=<P'(x 1 )J 1 (0(x 1 )) 



(j)'(x l ) ' 



Therefore, we can apply a transformation of the form (2.8) to arrange that Ji^x 1 ) = 0, and 
hence J = jo- Moreover, coordinates for which J is constant are uniquely determined up to an 
affine transformation 

x 1 = ax 1 + b, x 2 = ax 2 . 

Now suppose that a metric on the point-affine structure is given by 

Q (t* + uv 2 ) = Q [x 2 (A + ,o^) + « ) ) = \g(x)u 2 . (2.9) 

This case differs from the previous case in that the control vector field v% is already canonically 
defined by the point-affine structure prior to the introduction of a metric. Therefore, in order 
that the metric (2.9) be homogeneous, the unit control vector field 

1 

V2 



must be a constant scalar multiple of V2- Thus we must have G(x) = go for some positive 
constant go, and the homogeneous metrics in this case are given by quadratic functionals of the 
form 

Q(vi + uv 2 ) = ^9ou 2 
for some positive constant go, where v\, v 2 are the canonical frame vectors (2.4). 



2.2 Three states and one control 

In [1], we found three nontrivial local normal forms under point-affine equivalence. 

Remark 2.1. This classification assumes that the point-affine distribution is either bracket- 
generating or almost bracket-generating; otherwise the 3-manifold M can locally be foliated by 
a 1-parameter family of 2-dimensional submanifolds such that every trajectory of 3" is contained 
in a single leaf of the foliation. 

Case 2.1. 3~ = (gfr + x 3 ^ + J^) + span (^), The framing 

d 3 d ^ d d d j d 

1 dx l X dx 2 <9x 3 ' V2 dx 3 ' V3 Vl ' V2 q x 2 x3 dx 3 

(well-defined up to dilation in the (v 2 , i^-plane) has dual coframing 

r] 1 = dx 1 , rj 2 = dx 3 - J dx 1 - J x s (dx 2 - x 3 dx 1 ) , rf = dx 2 - x 3 dx 1 , 
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with structure equations 
drj 1 = 0, 

drj 2 = Ti^i] 1 A rf mod r] 2 , 
drj 3 = r/ 1 A rj 2 mod rf . 

As in Case 1.1, the method of equivalence does not lead to a completely determined coframing, 
so it is not clear from these structure equations whether this example is homogeneous as a point- 
affine distribution. 

So, suppose that a metric on the point-affine structure is given by 

Q ({w + + J i*) + « 5?) " 5 G M» 2 - < 2 ' 10 > 

The addition of the metric (2.10) allows us to choose a canonical framing (up to sign) by requi- 
ring V2 to be a unit vector for the metric, i.e., 

- 1 9 
and setting 

^3 = -[V1,V 2 ]. 

The canonical coframing associated to this framing is given by 

7 1 1 = dx 1 , rj 2 = ^G(x)(dx 3 - Jdx 1 ) modrj 3 , if = y/G(x) {dx 2 - x 3 dx 1 ) . (2.11) 

In order to identify the homogeneous examples, we consider the structure equations for the 
coframing (2.11), taking into account the fact that local coordinates for which the coframing 
takes the form (2.11) are determined only up to transformations of the form 

x 1 = x 1 +a, x 2 = 4>{x l ,x 2 ), x 3 = ^(x 1 ,^ 2 ) +x 3 </> i2 (x 1 ,x 2 ), (2.12) 

with ^ 0. Under such a transformation we have 



'G(x 1 ,x 2 ,5?') = ^(xV.a 3 )^, (2.13) 
J (x 1 , x 2 , x 3 ) = — ( J(x\ x 2 , x 3 ) - fa 2£ 2 (x 3 ) 2 - 2(j)^2X 3 - fei) , (2.14) 

with x , x 2 , x 3 as in (2.12). 

First consider the structure equation for drj 3 . A computation shows that 

7 3 ^x^ 2 3 i 1 

V = 2CP/ 2 ~ T1 A V m ° V ' 

Q „ 

Therefore, homogeneity implies that 2G x 3 / 2 must be equal to a constant —c\. The remaining 
analysis varies considerably depending on whether c\ is zero or nonzero. 
Case 2.1.1. First suppose that c\ = 0. Then G x 3 = 0, and so 

G(x\x 2 ,x 3 ) = G {x\x 2 ) 

for some function Go(x 1 ,x 2 ). According to (2.13), by a local change of coordinates of the 
form (2.12) with <j) a solution of the PDE 

M* 1 .^) = Go (fl, 0(^1^2))' 
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we can arrange that Gq^x 1 , x 2 ) = 1. This condition is preserved by transformations of the 
form (2.12) with 

<p(x l ,x 2 ) = x 2 + <j) {x 1 ). (2.15) 
With the assumption that G( X j X j X 3 ) = 1, the equation for drf reduces to 

drf = 7] 1 Ar] 2 + J^r] 1 A if. 
Therefore, J x 3 must be equal to a constant C3, and so 

J(x 1 ,x 2 ,x 3 ) = c 3 x 3 + Jq(x 1 , x 2 ) 
for some function Jq(x , x 2 ). Now the equation for drf becomes 

drf = (J ) X 2 rf A if . 
Therefore, (Jq) x 2 must be equal to a constant c 2 , and so 

J (x 1 , x 2 ) = c 2 x 2 + J\ (x 1 ) 
for some function J^x 1 ). With <\> as in (2.15) and 

J(x 1 ,x 2 ,x 3 ) = c 2 x 2 + c 3 x 3 + Jx^x 1 ), 
equation (2.14) reduces to 

J x {x l ) = Ji{x l +a) - (^'(z 1 ) - cs^x 1 ) - c 2 <f>o[x 1 )) . 

Therefore, we can choose local coordinates to arrange that Ji (x 1 ) = 0. 
To summarize, we have constructed local coordinates for which 

G(x\x 2 ,x 3 ) = 1, J(x 1 ,x 2 ,x 3 ) = c 2 x 2 + c 3 x 3 . 
These coordinates are determined up to transformations of the form 

x 1 = x 1 + a, x 2 = x 2 + (f>o (x 1 ) , x 3 = x 3 + (f)' (x 1 ) , 
where ^(x 1 ) is a solution of the ODE 

^'(x 1 ) - c^'^x 1 ) - c 2 <j)o(x 1 ) = 0. 
Case 2.1.2. Now suppose that c\ 7^ 0. Then 

G^x ,x , x 3 ^ = 2 

(cix 3 + Go(x 1 , x 2 )) 

for some function Gofx^x 2 ). According to (2.13), by a local change of coordinates of the 
form (2.12) with <j) a solution of the PDE 

4,1 (x\x 2 ) = —G (xS^x 1 ,^ 2 )) , 

we can arrange that Gq(x ,x ) = 0. This condition is preserved by transformations of the 
form (2.12) with 

<P(x\x 2 ) =<p (x 2 ). (2.16) 
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With the assumption that (^(x 1 , x 2 , x 3 ) = j—^s\s, the equation for drj 3 reduces to 

o ~t o (2</ x 3 J ^3) 19 99 

drf = 77 A 77 — 5 — — r? Art — c\jf A 7/ d . 

x- 1 

Therefore, ^g^ 3 ^ must be equal to a constant C3, and so 

J(x 1 ,x 2 ,x 3 ) = c 3 x 3 + J (x 1 ,x 2 )(x 3 ) 2 
for some function ^(x^x 2 ). Now the equation for drj 2 becomes 

drj 2 = -x 3 (J ) x i r, 1 At] 3 . 
The quantity — x 3 (Jo) x i can only be constant if (Jo) x i = 0; therefore, we must have 

J (x\x 2 ) = Ji(x 2 ) 
for some function Ji(x 2 ). With (j) as in (2.16) and 

J(x\x 2 ,x 3 ) =c 3 x 3 + Ji(x 2 )(x 3 ) 2 , 
equation (2.14) reduces to 

-„2\ Tlx (z.2\\ il (~2\ ^o(^ 2 ) 



0' o (x 2 )- 

Therefore, we can choose local coordinates to arrange that Ji(x 2 ) = 0. 
To summarize, we have constructed local coordinates for which 

G(x\x 2 ,x 3 ) = — 2, J(x 1 ,x 2 ,x 3 ) = c 3 x 3 . 

(cix 3 ) 

These coordinates are determined up to transformations of the form 

11 Q Q 

x = x + a, x = bx + c, x = bx + a 

Case 2.2. 3" = (x 2 g|x + x 3 ^ + J (x 2 ^)) + span (gfs-). We found a canonical framing 

2 d 3 <9 ./ , 3 \ 
ax 1 ox z \ ox 6 J 



2 d 

OX 6 



v 3 = ~[vi,v 2 ] = x 2 ^ + ((x 2 ) 2 J x3 - x 3 ) A (2.17) 
with dual coframing 



r, 1 - , 
x z 



1,9 1 T ,1 / T X 3 \/, 9 X 3 

(x 2 )^ 



X 2 X 2 \ X (t2\ 2 j \ X 2 



n" = -^.dx 3 ^Jdx 1 — I n I f dx 2 ^dx 1 

x z x A 

„3 _ 1 ,J2 _ & 

(x 2 )' 
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and structure equations 
dr] 1 = rj 1 A 7] 3 , 

drf = r] 1 A rj 2 + Tf z r] 1 A rf . (2.19) 

The local coordinates in the coframing (2.18) are only determined up to transformations of 
the form 

x 1 = cf)(x 1 ), x 2 = (j)'(x 1 )x 2 , x 3 = <p'(x 1 )x 3 + <P"(x r )(x 2 ) 2 , (2.20) 
with <p' (x 1 ) 7^ 0. Under such a transformation we have 

J(x\x 2 ,x 3 ) = J{x\x\x 3 ) - ^ (r^Wf + 3</>"(x 1 )x 3 ) , (2.21) 

with x 1 , x 2 , x 3 as in (2.20). 

First consider the structure equation for r/ 3 . Substituting the expressions (2.18) into the 
structure equation (2.19) for drj 3 shows that 

rp2 _ 2 t q^ 3 
J 12 — x J x 3 — o—^. 

x z 

Homogeneity implies that T 2 2 must be equal to a constant o, from which it follows that 

J(x\x 2 ,x 3 ) = | +a^ + J (x\x 2 ) 

for some function Jq(x ,x ). Now the equation for drj 2 yields 

T i3 = x 2 (Jo)x* - 2J - , 

and homogeneity implies that T 2 3 must be constant. The quantity (x 2 {Jo) x 2 — 2Jq — afjj) can 
only be constant if a = 0; therefore, we must have a = and 

x 2 {Jq) x 2 - 2J = -2ci 
for some constant c\. Therefore, 

J {x\x 2 ) =c 1 + J 1 (x 1 )(x 2 ) 2 
for some function Jjfx 1 ), and 



w 1 2 3\ ^ ( ^ ^ 



) 2 + c 1 + J 1 (x 1 )(x 2 ) 2 . 



With (j) as in (2.20) and J as above, equation (2.21) reduces to 
- / in // ix2 , , i \ \ (^'"(x 1 ) 3 ^"(x 1 ) 
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Therefore, we can choose local coordinates to arrange that J\ (x 1 ) = 0. This condition is pre- 
served by transformations of the form (2.20) with 



This implies that cp is a linear fractional transformation, i.e., 
v ' cx 1 +d 

Now suppose that a metric on the point-affine structure is given by 

Q(v 1 + uv 2 ) = \g{x)u 2 . (2.22) 

As in Case 1.2, the control vector field v 2 is already canonically defined by the point-affine 
structure prior to the introduction of a metric. Therefore, in order that the metric (2.22) be 
homogeneous, the unit control vector field 

1 

: V 2 



y/G(p 

must be a constant scalar multiple of v 2 . Thus we must have G{x) = go for some positive cons- 
tant go, and the homogeneous metrics in this case are given by quadratic functionals of the form 

Q(vi + uv 2 ) = l^gou 2 

for some positive constant go, where vi, v 2 , V3 are the canonical frame vectors (2.17). 
To summarize, we have constructed local coordinates for which 



,3\ 2 



G(x\x 2 ,x 3 ) = go, J(x\x 2 ,x 3 ) = I fjj^l +ci. 

These coordinates are determined up to transformations of the form 

1 ax 1 + b 9 ad — be 9 q ad — be o 2c(ad — be) 9 

r^ 1 — rg^ — — V — X 

cxi + cf (c^+d) 2 ' (cx x + d) 2 (c^+d) 3 
Case 2.3. 

J - (— j(x 3 — 9 I g 9 \] +s an fx 3 — — + H — 

ydx 1 \ Sx 1 <9x 2 <9x 3 / / S ^ \ Sx 1 <9x 2 <9x 3 

where 7^ 0. We found a canonical framing 

v - 9 + J fx 3 9 + 9 + H 9 
dx 1 \ dx 1 dx 2 dx c 

e ( x zA_ , A + # JM 
y/eH x i \ dx 1 dx 2 dx 3 J ' 



^3 = -[Ul, «2] 
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where e = ±1 = sgn(H x i), with dual coframing 
r/ 1 = dx 1 — x 3 dx 2 , 



7] 2 = e^/eH x i (dx 2 - j(dx 1 - x 3 dx 2 )) mod rf ', 

rf = - 1 (H dx 2 - dx 3 ) , (2.23) 

and structure equations 

dr, 1 =Tltf X Ar, 3 +T^ri 2 Arf, 
dr ] 2 =T 2 3 r ] 1 Ar ] 3 +T 2 3 r ] 2 Ar ] 3 , 

dr] 3 = v 1 Ar] 2 + T^r] 1 A rf + T^rf A rf. (2.24) 

The identification of homogeneous examples is considerably more complicated than in the pre- 
vious cases. We refer the reader to Appendix A for the details. We find that the homogeneous 
examples in this case are all locally equivalent to one of the following: 

• J(x 1 , x 2 , x 3 ) = cx, H{x x ,x 2 , x 3 ) = e(x l + C2X 3 ) 
for some constants c\, C2; 

• J(x 1 ,x 2 ,x 3 ) = a cos (c 3 x l )/\J ec 3 (c 3 (x 3 ) 2 + c 4 ), 

#(x\:r 2 ,:r 3 ) = (c 3 (x 3 ) 2 + c 4 ) tan (c^x 1 ) + F 20 (x 2 ) J c 3 (x 3 ) 2 + c 4 

for some constants c\, c 3 , c 4 with c 3 ^ 0, and some arbitrary function F2o{x 2 )] 

• J(x 1 ,x 2 ,x 3 ) = a cosh (C3X 1 ) I ' \j ec 3 (c 3 (x 3 ) 2 -c 4 ), 

H(x l ,x 2 ,x 3 ) = (-c 3 (x 3 ) 2 + c 4 ) tanh (c 3 x l ) + F 2 o(x 2 ) \J c 3 (x 3 ) 2 - c 4 

for some constants ci, c 3 , c 4 with c 3 ^ 0, and some arbitrary function F2o(x 2 ). 

Now suppose that a metric on the point-affine structure is given by 
Q (vi + UV2) = ^G(x)u 2 . 

As in the previous case, since the control vector field V2 is already canonically defined by the 
point-affine structure prior to the introduction of a metric, we must have G(x) = go for some 
positive constant go- 

The results of this section are encapsulated in the following two theorems: 

Theorem 2.2. Let 3" be a rank 1 strictly affine point-affine distribution of constant type on 
a 2- dimensional manifold M , equipped with a positive definite quadratic cost functional Q. If 
the structure (3~, Q) is homogeneous, then (3~, Q) is locally point-affine equivalent to 

J = vx + span (v 2 ), Q{vx + uv 2 ) = ^G(x)u 2 , 
where the triple (vx,V2,G(x)) is one of the following: 
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Theorem 2.3. Let J be a rank 1, strictly affine, bracket-generating or almost bracket- generating 
point-affine distribution of constant type on a 3 -dimensional manifold M, equipped with a positive 
definite quadratic cost functional Q . If the structure (3~, Q) is homogeneous, then (3~, Q) is locally 
point-affine equivalent to 

3" = v\ + span (v 2 ), Q(vi + uv 2 ) = ^G(x)u 2 , 
where the triple (vi,v 2 ,G(x)) is one of the following: 

(2- L1 ) ^ = Zi + xZ 4j2 + { c 2X 2 + c^)^ "2 = ^, G(x) = l; 



dx 1 dx 2 dx 3 ' dx 3 

(2.1.2) = i + ^ 3 i +c 3 x 3 ^, v 2 = ^, G{x 
OX L ox z ox° ox a 



(cix 3 ) 



2 ' 



x 2 J J \ dx 3 



v 2 = x 2 -^, G{x) = g ; 
(2-3-1) v x = A + Cl (x^ + ^ + e(x^ + c 2 * 3 ) ^ ) , 



v 2 = e ( x 3 A + A + e(x i + C2X 3 } ^ G{x) = go . 



d ci cos(c3X 1 ) / 3 d d d \ 

{ ] Vl = dx^ + 7^3(^)2 + C4) [ X d^ + d^ +H dx- 3 )> 

U2 = T^ + ^ + i ^J' G ^ = 9 ^ 

where H = ^(c 3 (x 3 ) 2 + c 4 ) tan (C32; 1 ) + F 20 (x 2 ) ^ c 3 (x 3 ) 2 + c 4 ^ ; 
9 cicosh(c3X 1 ) / 3 <9 3 3 \ 

C2-3-3) ^ = ^t + 7^3^3)2 _ C4) V x a^ + ^ + // ax3j' 

U2 = T^ + ^ + i ^J' G ^ = 50 ' 

w/iere il = ^( — c 3 (x 3 ) 2 + c 4 ) tanh (C3X 1 ) + -F20 (^ 2 ) \J C3 (z 3 ) 2 - c^j . 

3 Optimal control problem for homogeneous metrics 
3.1 Two states and one control 

In this section we use Pontryagin's maximum principle to compute optimal trajectories for each 
of the homogeneous metrics of Theorem 2.2. 

Case 1.1. This point-affine distribution corresponds to the control system 

x l = 1, x 2 = u, (3.1) 

with cost functional 

Q(x) = \e 2 ^\ 2 . 
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Consider the problem of computing optimal trajectories for (3.1). The Hamiltonian for the 
energy functional (1.2) is 

U =pi± l + p 2 x 2 -Q(x) =pi +p 2 u- ^e 2cixl u 2 . 



By Pontryagin's maximum principle, a necessary condition for optimal trajectories is that the 
control function u(t) is chosen so as to maximize %. Since u is unrestricted and ^e 2ciX > 0, 
m&x u l-L occurs when 



= ^=p 2 -e 2 ^u, 

that is, when 

u = p 2 e- 2c ^\ 
So along an optimal trajectory, we have 

H =Pi + {P2?e~ 2 ^ - \{p 2 ?e- 2 ^ = P1 + \{p 2 ) 2 e- 2 ^\ 
Moreover, Ti is constant along trajectories, and so we have 
Pi + \{p2?e- 2c ^ =k. 

Hamilton's equations 

. _&H . _ _&H 

dp ' ^ dx 



take the form 



x l = l, p 1 = Cl ( P2 ) 2 e- 2c ^ , 

x 2 = p 2 e~ 2cixl , p 2 = 0. (3.2) 



The equation for p 2 in (3.2) implies that p 2 is constant; say, p 2 = c 2 . Then optimal trajectories 
are solutions of the system 

x 1 = l, x 2 = c 2 e- 2cix \ 

This system can be integrated explicitly: 

• If c\ = 0, then the solutions are 

x 1 = t, x 2 = c 2 t + C3. 

These solutions correspond to the family of curves 

x 2 = c 2 x x + c 3 

in the fx 1 , x 2 ) -plane. Thus, the set of critical curves consists of all non-vertical straight 
lines in the (x 1 ,^ 2 ) plane, oriented in the direction of increasing x 1 . 
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• If c\ 7^ 0, then the solutions are 

x x = t, x 2 = -— c 2 e- 2ci *. 

2ci 

These solutions correspond to the family of curves 
x 2 = ~C2e~^ 

in the (x 1 , x 2 ) -plane. Thus, the set of critical curves consists of a family of exponential 
curves in the fx 1 , a; 2 ) plane, oriented in the direction of increasing x 1 . 

Case 1.2. This point-affine distribution corresponds to the control system 

■ 1 2 -2 2 • i 2 

X = X , X = X Jo + X U, 

with cost functional 
Q( x ) = T^dou 2 - 

Pontryagin's maximum principle leads to the Hamiltonian 
U = pix 2 + p 2 x 2 j + 7— (p 2 x 2 ) 2 

2#0 

along an optimal trajectory, and Hamilton's equations take the form 
x 1 = x 2 , pi = 0, 

■2 2- 1 P2(x 2 ) 2 . . {P2) 2 x 2 

x =xj H — , P2 = -Pi-P2jo • (3.3) 

9o 9o 

It is straightforward to show that the three functions 

h = U = p\x 2 +p 2 x 2 j + ^—(p 2 x 2 ) 2 , h=Pl, I 3 =pix 1 +p 2 x 2 

are first integrals for this system. This observation alone would in principle allow us to construct 
unparametrized solution curves for the system. But in fact, we can solve this system fully, as 
follows. 

The equation for p\ in (3.3) implies that p\ is constant; say, p\ = c\. Now it is straightforward 
to show that 

j f {p 2 x 2 ) + Cl x 2 = 0. (3.4) 

If ci = 0, then (3.4) implies that p 2 x 2 is equal to a constant k 2 , and so 

■ 2 2 / ■ , k 2 \ 2 
x =x jo H ) =c 2 x . 

V 90 J 

There are two subcases, depending on the value of c 2 . 

• If c 2 = 0, then x 2 = C3, and since x 1 = x 2 , we have x 1 = c%t + C4. These solutions 
correspond to the family of curves x 2 = C3 in the (x 1 , x 2 ) -plane. These curves are all 
horizontal lines, oriented in the direction of increasing x 1 when x 2 > and decreasing x 1 
when x 2 < 0. 
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• If C2 7^ 0, then x 2 = c 3 e C2 *, and since x 1 = x 2 , we have x 1 = j|e C2 * + c 4 . These solutions 
correspond to the family of curves x 2 = C2(x 1 — cA in the (x 1 , x 2 ) -plane. These curves 
are all non-vertical, non-horizontal lines, oriented in the direction of increasing x 1 when 
x 2 > and decreasing x 1 when x 2 < 0. 

On the other hand, if c\ 7^ 0, then it is straightforward to show that 

d 2 J , 



I 2^ d , 2\f ■ , P2X 2 \ 



dt 2 

Integrating this equation once gives 

|(^ 2 )=Jo(p 2 x 2 ) + ^gS + c 2 . (3.5) 

There are three subcases, depending on the value of k = go(jo9o — 2C2). 

• If k = 0, then the solution to (3.5) is 

2 go(2 + io(t + c 3 )) 

P2X = — , 

t + c 3 

and from equation (3.4), 

(p 2 x ) ~ 



c\ dt ' ci(t + c 3 ) 2 ' 

Then since x 1 = x 2 = — ^-^(P2% 2 ), we have 

1 1 / 2\ , 5o(2 + jo(t + c 3 )) 
x = (p 2 x z ) + c 4 = — r hc 4 . 

These solutions correspond to the family of curves 



x 



2 



(cix 1 - (j g + C1C4)) 2 



2cigo 



in the (x 1 , x 2 ) -plane. These curves are all parabolas with vertex lying on the a^-axis. 
Since we must have x 2 7^ 0, the set of critical curves consists of all branches of parabolas 
with vertex on the ir 2 -axis, oriented in the direction of increasing x 1 when x 2 > and 
decreasing x 1 when x 2 < 0. 

If k > 0, then the solution to (3.5) is 

p 2 x 2 = -V&tanh (^~( t + c ^ ~hgo, 
and from equation (3.4), 

„2 _ 1 d , 2 \ _ k „ n 2 ' 



(p 2 x ) = sech - — (t + c 3 ) 



cidt K 2cig y2g 

Then since x 1 = x 2 = — ^-^(P2X 2 ), we have 
.1 _ 1 f„ ~2\ , „ _ 1 L/r. , / 



' ~ - — (P2^c 2 ) + c 4 = — I \/fctanh ( — (t + c 3 ) j + jo5o j + c 4 . 
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These solutions correspond to the family of curves 
x 2 = [{cix 1 - (j g Q + cic 4 )) 2 - k] 

in the (x 1 , x 2 ) -plane. These curves are all parabolas opening towards the x 1 -axis. Thus the 
set of critical curves consists of parabolic arcs opening towards the x 1 -axis, approaching 
the axis as t — > ±00, and oriented in the direction of increasing x 1 when x 2 > and 
decreasing x 1 when x 2 < 0. 

If k < 0, then the solution to (3.5) is 

p 2 x 2 = V^ktan (^^-(t + c 3^ -Jo9o, 
and from equation (3.4), 



Then since x 1 = x 2 = ~^^{P2X 2 ) , we have 



x 

c\ dt ' 



x 1 = -^-(p 2 x 2 ) +c 4 = ^v^tan (~^~~( t + c 3)^ ~ 3o9o) + c 4 . 
These solutions correspond to the family of curves 
x 2 = — t~ \{c\x l - (j g + cic 4 )) 2 - k] 

in the (x , x 2 )-plane. These curves are all parabolas opening away from the x 1 -axis. Thus 
the set of critical curves consists of parabolic arcs opening away from the x 1 -axis, becoming 
unbounded in finite time, and oriented in the direction of increasing x 1 when x 2 > and 
decreasing x 1 when x 2 < 0. 

3.2 Three states and one control 

In this section we use Pontryagin's maximum principle to compute optimal trajectories for each 
of the homogeneous metrics of Theorem 2.3. 

Case 2.1.1. This point-affine distribution corresponds to the control system 

X = 1, X =X , X = C2X + C3X + u, 

with cost functional 
Q(x) = l -u 2 . 

The Hamiltonian for the energy functional (1.2) is 

% = Pl± +p 2 X +P3X -Q(x)=pi+P2X +P3(c 2 X +C3X +U)--U. 

Pontryagin's maximum principle leads to the Hamiltonian 
U=pi+ p 2 x 3 + p 3 (c 2 x 2 + c 3 x 3 ) + - (p 3 ) 2 
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Figure 1. 

along an optimal trajectory, and Hamilton's equations take the form 

x 1 = 1, p\ = 0, 

x 2 = x 3 , p 2 = ~C 2 P3, 

x 3 = c 2 x 2 + c 3 x 3 +p 3 , p 3 = -p 2 - c 3 p 3 . (3.6) 
The equations for p 2 and p 3 in (3.6) can be written as 
P2 + c 3 p 2 - c 2 p 2 = 0, 

and the function p 3 = —^P2 satisfies this same ODE. Then the equations for x 2 and x 3 can be 
written as 

x -c 3 x -c 2 x =p 3 {t), 
where p 3 (t) is an arbitrary solution of the ODE 

P3 + c 3 p 3 - c 2 p 3 = 0. 
Therefore, x 2 (t) is an arbitrary solution of the 4th-order ODE 

/ d 2 d \ ( d 2 d \ 2 . . 

{d^ +C3 dt- C2 ) {dri- C3 dt- C2 ) x{t) = > 

and for any such x 2 (t), we have 

x 1 (t)=t + t , x 3 (t)=x 2 (t). 

A sample optimal trajectory is shown in Fig. 1. 

Case 2.1.2. This point-affine distribution corresponds to the control system 

X = 1, X =X , X = C 3 X +u, 

with cost functional 

Q{±) = 1 u 2 . 
2[cix 6 ) 
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Figure 2. 



Figure 3. 



Pontryagin's maximum principle leads to the Hamiltonian 
U=pi + P2X 3 + c 3 p 3 x 3 + - (ax 3 p 3 ) 2 
along an optimal trajectory, and Hamilton's equations take the form 



x l = 1, 



Pi = 0, 

C 3 = C 3 X 3 + (ciX 3 ) 2 p3, p 3 = ~P2 ~ C 3 p 3 - (c\p^ 2 



X 2 = X 3 



X . 



(3.7) 



The equation for p2 in (3.7) implies that pi.it) is equal to a constant C2- Then (3.7) implies that 
(p 3 x 3 ) = -c 2 x 3 , x 3 = c 3 x 3 + cix 3 (p 3 x 3 ). (3.8) 
These equations can be solved as follows: 

• If C2 = 0, then the function p 3 x 3 is constant, and so the equation for x 3 becomes 



x 3 = cx 3 



for some constant c. If c = 0, then the solution trajectories are given by 

x 1 (t) = t + t , x 2 (t) = at + b, x 3 (t) = a 

for some constants a, b. Sample optimal trajectories are shown in Fig. 2. 
If c 0, then the solution trajectories are given by 

x l {t) = t + t Ql x 2 (t) = a e a + b, x 3 (t) = ae dt 

c 

for some constants a, b. Sample optimal trajectories are shown in Fig. 3. 

If C2 7^ 0, then (3.8) can be written as the 2nd-order ODE for the function z(t) = p 3 (t)x 3 (t): 



(c 3 + c\z) 



z. 
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Figure 4. 

Integrating once yields 

z = -^(ciz) 2 + C 3 Z + c 4 

for some constant C4. Depending on the values of the constants, the solution z(t) has one 
of the following forms: 



(1) z (t) = a tan(6t + c) + d, 

(2) z{t) = a tanh(5t + c) + d, 

(3) zCH^ + c, 



c 3 — 2C1C4 < 0; 
c 3 — 2C1C4 > 0; 

Co — 2C1C4 = 0. 



Then we have 



-z = X 



C2 



and so the corresponding solution trajectories are given (with slightly modified constants) 
by: 



(3) 



'x l {t) = t + t , 
x 2 (t) = atan(6t + c) + d, 
^x 3 (t) = ab sec 2 (bt + c); 

= t + t , 

x 2 (t) = — — + c, 

x 3 (t) 



at + 6 



(at + 6) 



2 ' 



(2) 



V(t) = t + t , 
x 2 (t) = atanh(6t + c) + d, 
^x 3 (t) = a6sech 2 (6t + c); 



Sample optimal trajectories for the first two cases are shown in Fig. 4. 

Case 2.2. This point-affine distribution corresponds to the control system 



x 1 = x 2 , 



x 2 = x 3 , 



x 3 = x 2 



3 fx 



3\ 2 



2 \x 



+ Cl+U 
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with cost functional 



Q( x ) = 7,90V 2 . 



Pontryagin's maximum principle leads to the Hamiltonian 

2 s 2 l 3 / x 3 \ 2 \ 1 , 2 s 2 

ri = pix + p 2 x + p 3 x I - I +ci j + ^{P3X ) 

along an optimal trajectory, and Hamilton's equations take the form 



x 1 = x 2 , pi = 0, 



■2 3 • <j yjy-^ j x / \2 2 

X =X, p 2 = ~Pl + ~ 2 ClP3--(P3)X 



3p 3 (> 3 ) 2 1 



X 3 



3 (x 3 ) 2 2 1 / 2 \2 P3X 3 

2~^ + ClX + g P3 ( x ) ' ^3 = -P2-3-^2-. (3.9) 



The system (3.9) has three independent first integrals in addition to the Hamiltonian Ti 
(which is automatically a first integral): it is straightforward to show, using (3.9), that the 
functions 

h=Pl, h = Pix 1 +P2X 2 + p 3 x 3 , I 3 = pi(x 1 ) 2 + 2p 2 x 1 x 2 + 2p 3 x 1 x 3 + 2p 3 (x 2 ) 2 

are first integrals for this system. We can use these conserved quantities to reduce the sys- 
tem (3.9), as follows: on any solution curve of (3.9), we have 

h = h, h = k 2 , h = k 3 

for some constants k%, k 2 , k 3 . These equations can be solved for p\, p 2 , p 3 to obtain 

Pi = h, 

, / x 1 (x 1 ) 2 ^ 3 ^ , / 1 x 1 x 3 \ , / x 3 
P2 = h --5 - -. „.„ + k 2 -5 + )+k 3 ' 



x 2 2(x 2 ) 3 J "\x 2 (x 2 ) 3 / V 2 (^ 2 ) 3 / ' 



P3 = ki[ 777-^ + k 2 [ - j-^ + k. 



2(x 2 ) 2 7 V (^ 2 ) 2 J \2(x 2 ) 2 

These equations can be substituted into (3.9) to obtain a closed, first-order ODE system for 
the functions X . X . X , depending on the parameters k\, k 2 , k 3 ; moreover, making the same 
substitution in the Hamiltonian Ti yields a conserved quantity for this system. (The precise 
expressions for the system and the conserved quantity are complicated and unenlightening, 
so we will not write them out explicitly here.) The resulting ODE system cannot be solved 
analytically, but numerical integration yields sample trajectories as shown in Fig. 5. 
Case 2.3.1. This point-affine distribution corresponds to the control system 

x 1 = 1 + x 3 (ci + u), x 2 = c\ + u, x 3 = e(x 1 + C2X 3 ) (c\ + u), 

with cost functional 

Q(x) = ^gov 2 - 
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Figure 5. 



Pontryagin's maximum principle leads to the Hamiltonian 

W=Pl+ P3X 1 + -7^ (piX 3 +p 2 + C 3 p 3 X 3 ) 
V c 2 



+ (piX 3 +P2 + C 3 p 3 X 3 ) (c 2 (piX 3 +p 2 +P3X 1 + C 3 p 3 X 3 )) 

along an optimal trajectory, and Hamilton's equations take the form 
x 1 = 1 + x 3 ( + c 2 (pix 3 +P2+ P3X 1 + c 3 p 3 x 3 ^ 

■ 2 c l 

x — 



+ C 2 (piX 3 + p 2 + + C 3 p 3 X 3 ) , 



X 3 = (x 



+ C 3 X 3 ) + C 2 (piX 3 +p 2 +p 3 X l + C 3 p 3 X 3 )^ , 

Pi = ~P3 +C 2 (piX 3 +p 2 +P3.X 1 +C 3 p 3 X 3 ) S j , 



P2 = 0, 

P3 = ~ ^(Pl + C3P3) (^^= + C 2{P\X 3 +P2 +P3X 1 + C 3 p 3 X 3 )^ 

+ {(Pi + c 3 p 3 )x 3 + p 2 ) [C2P\ + c 2 c 3 p 3 )^ . (3.10) 

The system (3.10) has three independent first integrals in addition to the Hamiltonian Ti: it 
is straightforward to show, using (3.10), that the functions 

h = (pi + np 3 )e rix2 , I 2 = (pi + r 2 p3)e r2X \ I 3 = p 2 , 

where r\, r 2 are as in (A. 19), are first integrals for this system. A similar process to that 
described in the previous case leads to a closed, first-order system of ODEs for the func- 
tions x , x 2 , x 3 ; numerical integration of this system yields sample trajectories (with c 3 > 
and c 3 < 0) as shown in Fig. 6. 

Case 2.3.2. Due to the complexity of the computations, we will restrict our attention to the 
simplest case, where e = 1, c\ = C4 = 0, and -F 2 o(^ 2 ) = 0. (In the interest of brevity, we will 
omit Case 2.3.3, which is similar to this case.) With these assumptions, we have 

d n d d / 3 \2 / n 9 

vi = ^-r, v 2 = x —- r + —~ + c 3 (x J tan (c 3 x l-^-s-. 
ox L ox 1 ox z ' ' ox 6 
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Figure 6. 



This point-affine distribution corresponds to the control system 



X = 1 + X U, X = u, 



with cost functional 



c 3 (x 3 ) 2 tan {c^x l ^ju, 



Q( x ) = ^9ou 2 - 

Pontryagin's maximum principle leads to the Hamiltonian 



1 '2ff^A*^2f^ , o-.„ -Wf^n+pfitnLfax 1 ) 



H = W~ (4{{x 3 ) Ps) tan 2 (csx 1 ) + 2c 3 p 3 
^9o v 

+ (2 P1 (x 3 ^ + So) + {x 3 Pi) 2 + vl 



along an optimal trajectory, and Hamilton's equations take the form 

J ~ ~o ( C3 ( x3 ) 3 ^ 3tan ( c 3 xl ) + ( x3 ) 2 pi + x3 P2 + 9o) , 



x 



x 2 = — (c 3 (x 3 ) 2 p 3 tan (C3X 1 )) + x 3 p 1 +p 2 , 



x 3 



Pi 



1 c 3 (x 3 ) 2 tan (C3X 1 ) (c 3 (x 3 ) 2 p 3 tan (c^x 1 ) + x 3 pi + p 2 ) , 



90 



9o cos° 



7 ry (c 3 x 3 ) 2 p 3 (c 3 (x 3 ) 2 p 3 sin (C31 1 ) + (x 3 pi + p 2 ) cos (c 3 x l 

( i C 3 x i J v 



P2 = 0, 



1 



P3 = -— (2(c 3 p 3 ) 2 (x 3 ) tan 2 (qx 1 ) + c 3 x 3 p 3 (3x 3 pi + 2p 2 ) tan (qx 1 ) +pi(x 3 pi + p 2 

This ODE system cannot be solved analytically, but numerical integration yields sample trajec- 
tories as shown in Fig. 7. 



4 Conclusion 



What is perhaps most interesting about these results is how the behavior of control-affine 
systems in low dimensions varies from that of control-linear (i.e., driftless) systems. As we 
observed in [1], functional invariants appear in much lower dimension for affine distributions 
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4' 



10 



Figure 7. 



(beginning with n = 2, s = 1) than for linear distributions, where there are no functional 
invariants in dimensions below n = 5, s = 2. 

With the addition of a quadratic cost functional, we see a similar phenomenon: for linear 
distributions with a quadratic cost functional, there are no functional invariants for any n when 
s = 1, since local coordinates can always be chosen so that a unit vector field for the cost 
functional is represented by the vector field J^. But for affine distributions with s = 1, there 
are numerous functional invariants, and even the homogeneous examples exhibit a wide variety 
of behaviors for the optimal trajectories. 

A Normal forms for Case 2.3 

In this appendix, we carry out the analysis to identify examples of normal forms for homogeneous 
point-affine structures in Case 2.3. 

First consider local coordinate transformations which preserve the expressions (2.23). Let 
fx 1 , x 2 , x 3 ) and (x 1 , x 2 , x 3 ) be two local coordinate systems with respect to which the coframing 
{f] 1 ,rj 2 ,rj 3 ) takes the form (2.23). Then we must have 



r] 1 = dx 1 



— x' 



■ 3 dx 2 = dx 1 - x 3 dx 2 . 



(A.l) 



Taking the exterior derivative of (A.l) yields 



dr/ 1 = dx 2 A dx 3 = dx 2 A dx 3 . 



(A.2) 



In particular 



span(cte 2 , dx 3 ) = span(dx 2 , dx 3 ) . 



Therefore we must have 





(A.3) 




Equation (A.2) then implies that the functions <j>, tp 



satisfy the PDE 



= 1. 



(A.4) 
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Unfortunately, equation (A. 4) cannot be solved explicitly in terms of arbitrary functions 
of x 2 , x 3 . However, it can be solved implicitly with a slightly different setup. Instead of (A. 3), 
suppose that we define our coordinate transformation by 

x 2 = <p(x 2 ,x 3 ), x 3 = ip(x 2 ,x 3 ). 

Then equation (A. 2) is equivalent to the condition 

(In addition, both terms in this equation must be nonzero.) This is equivalent to the condition 
that there exists a function $(x 2 ,x 3 ) such that 

0(X 2 ,X 3 ) = $53, ?P(X 2 ,X 3 ) = $ X 2. 

Then equation (A.l) implies that 

x 1 = x l + $(x 2 ,x 3 ) - x 3 $ x 3(x 2 ,x 3 ). 

The local coordinate transformations which preserve the expression for r/ 1 in (2.23) are defined 
implicitly by 

x 1 = x l + 3>(x 2 ,x 3 ) - x 3 $ x 3(x 2 ,x 3 ), x 2 = ^ x s(x 2 ,x 3 ), x 3 = $ x 2(x 2 ,x 3 ), (A. 5) 

where <I>(x 2 ,x 3 ) is an arbitrary smooth function of two variables with $ x 2%3 7^ 0. 

Next we will compute how the function if (x 1 , x 2 , x 3 ) transforms under a coordinate trans- 
formation of the form (A. 5). (When we consider the implications of homogeneity, it will turn 
out that J can be expressed in terms of H and its derivatives; thus there is no need to explicitly 
compute the effects of the transformation (A. 5) on J.) Consider the expression for rj 3 in (2.23). 
We must have 

if = 1 (H(x) dx 2 - dx 3 ) = 1 (H(x) dx 2 - dx 3 ) . (A.6) 
VWO*) \/eH^(x) 

From (A. 5), we have 

dx 2 = & X 2 X 3 dx 2 + & X 3 X 3 dx 3 , dx 3 = <& X 2 X 2 dx 2 + ® x 2 x 3 dx 3 . 

Substituting these expressions into (A.6) yields 

—====((H(x) - <$> X 2 X 2) dx 2 - <S> X 2 X 3 dx 3 ) 
^eH x i(x) 

1 (H(x)$ x 2 x 3 dx 2 + (H(x)$ x s x 3 - l)dx 3 ). (A.7) 



<//,:!•'•) 

Equating the ratios of the coefficients of dx 2 and dx 3 on both sides of (A.7) yields 

(H(x) -$ X 2 X 2) H{X)<S> X 2 X 3 



-*x»5* (F(x)cD^3-l) 
which implies that 



h(x\x z , x 3 ) = vv x x ' x x xx y (A 8) 

1 7 l-^H(x\x 2 ,x 3 ) 
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Now suppose that the structure is homogeneous. Unlike the previous cases, the assumption 
of homogeneity will imply some relations among the constants appearing in the structure equa- 
tions (2.24). In the homogeneous case, the functions Tj k are all constant, and differentiating 
equations (2.24) implies that 

= d{dr ] 1 ) = (T^Tfa - T 1 ^) r, 1 A rf A r? 3 , 
= d{drf) = (2*2* - 2f 3 2f 3 ) r, 1 A rf A r, 3 , 
= d(dr, 3 ) = -(2^ + T 2 3 ) ?/ A r, 2 A r, 3 . 

The first two equations imply that the vectors 

[ThT^l [Tf 3 T 2 3 ], [T 3 3 T 3 3 ] 

are all scalar multiples of each other unless T 3 3 = T| 3 = 0, while the third equation implies that 

1 23 — 1 13- 

In most of the computations that follow, these relations will be self-evident; however, at one 
point they will have implications for the function H. 
The structure equation for dr] 1 is 

dr/ 1 = -J^eH x ir] 1 A if - erj 2 A rf . 

Therefore, we must have 

J Cl 



for some constant c\, so that the equation for dr] 1 becomes 
dr] 1 = -ci?] 1 At/ 3 — erf A rf . 
Now the equation for drf reduces to 
drf = r] 1 Ar] 2 



' 'H x i x 2 + x 3 H x i x i + HH x i x3 - 2H x iH x z) (ar, 1 A rf + erf A rf) . 



22/3.1 \J eH x i 
Therefore, 

H x i x 2 -\- x 3 H x i x i + HH x i x s — 1H x iH 



3 



H x iy/eH x i 



-2c 2 (A.9) 



for some constant C2. Substituting the derivative of (A.9) with respect to x 1 into the equation 
for drf yields 

2 1 u . , . „2 



Observe that: 



The coefficient of rf A rf in drf is equal to minus the coefficient of rf A if in drf , as we 
previously observed that it must be. 
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• If C2 7^ 0, then the ratio of the rj 1 A rj 3 and rj 2 A r] 3 coefficients in drf must be equal to — 

(which is the ratio of these coefficients in dry 1 ), and hence the r/ 1 A rj 3 coefficient in drj 2 

(? 

must be equal to -f- 
Therefore, if c 2 ^ 0, then H satisfies the PDE 

The solutions of (A. 10) are precisely the linear fractional transformations in the x 1 variable, 
and so we must have 

, , 9 o N F 1 (x 2 ,x 3 )x l + F (x 2 ,x 3 ) 

H(x\x 2 ,x 3 ) = ) ' J. l , ' /. A.ll 

V ' ' 7 GiO 2 ,^ 1 + G (x 2 ,x 3 ) 1 7 

for some functions F (x 2 ,x 3 ), iq(x 2 ,x 3 ), Go(x 2 ,x 3 ), Gi(x 2 ,x 3 ). 

By contrast, if c 2 = 0, then the vectors [T 13 T 2 1 3 ], [T 2 3 T 2 2 3 ] are no longer required to be 
linearly independent, and so the Schwarzian derivative of H with respect to x 1 appearing in 
equation (A. 10) is only required to be constant, but not necessarily equal to zero. There are 
two possibilities, depending on the sign: 

• If § (%^) 2 - = -4 for c 3 > 0, then 

\ X 1 / X L 

H(x 1 ,x 2 ,x 3 ) = Fi(x 2 ,x 3 )tan(c 3 x 1 + F (x 2 ,x 3 )) + F 2 (x 2 ,x 3 ) (A.12) 
for some functions Fq[x 2 , x 3 ) , F\ (x 2 , x 3 ) , F2^x 2 ,x 3 ) with .Fi(x 2 ,x 3 ) ^ 0. 

• If f (%^) 2 - = c| for c 3 > 0, then 

\ X 1 / X L 

H(x 1 ,x 2 ,x 3 ) = Fi(x 2 ,x 3 )tanh(c 3 x 1 + F (x 2 ,x 3 )) +F 2 (x 2 ,x 3 ) (A.13) 
for some functions i*o(x 2 ,x 3 J, F\ (x 2 , x 3 ) , F2(x 2 ,x 3 ) with i*i(x 2 ,x 3 ) ^ 0. 
We consider each of these cases separately. 

A.l c 2 ^ 

In this case, Hfx 1 , x 2 , x 3 ) is given by (A.ll). Now we compute how the function (A.ll) trans- 
forms under a local coordinate transformation of the form (A. 5). 

Lemma A.l. There exists a local coordinate transformation of the form (A. 5) such that the 
function ^(x 1 ,^ 2 ,^ 3 ) is linear in x , i.e., 

H{x l ,x 2 ,x 3 ) = Fifx 2 ,^ 3 )^ 1 + F (x 2 ,x 3 ), (A.14) 
with Fi ^ 0. 

Proof. Equation (A. 8) can be written as 

H(x 1 ,x 2 ,x 3 ) = H(x 1 ,x 2 ,x 3 ) - $ X 2 X 2 ^ (A.15) 

®&&H{X 1 , X 2 , X 3 ) + ((<P x 2 x3 ) 2 - § x 2 x 2<S> x 3 x z) ' 
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Substituting (A. 11) into this equation yields 

~ , ( F 1 x 1 +F \ ^ 



Gi* 1 + Go $-3-3 (gfr±f) + ((fea) 2 - <IVv'<IV.,0 

^ (iV + Fp) - $ X 2 X 2 (dx 1 + Gp) 

**3x3 O^ 1 + F ) + (($ x 2 i3 ) 2 - $> x 2 x2 $ i3 z 3 ) (dx 1 + G ) 

= [Fl - § x 2 x 2Gl] X 1 + [Fp - $> X 2 X 2 Gp] 

" [^3£3Fl + (($ x 2 ; g3) 2 -$ :c 2 a; 2$£3 £ 3) G X ] X 1 + <l> /•'„ • i i <!>,,■,.. i 2 <l> «l> ) G ] ' 

The coefficients of x 1 on the left-hand side of this equation are the same as the coefficients of x 1 
on the right-hand side, so the condition that G\ = is equivalent to 



= ^^(x 2 ,^ 3 )^ 2 ^ 3 ) + (($ x 2 x3 (x 2 ,x 3 )) 2 - <5> x 2 x 2(x 2 ,x 3 )<S> x3x3 (x 2 ,x 3 )) Gi(x 2 ,x 3 ) 
= $ £3£3 (x 2 , x 3 )iq (x 2 , <$> X 2 (x 2 , x 3 ) ) 

+ ((^ 2£ 3(x 2 ,x 3 )) 2 -^2 x 2(x 2 ,x 3 )$^3(a; 2 ,x 3 )) Gi (x 2 , ® x 2 (x 2 , x 3 )) . 



Any solution <I>(x 2 ,x 3 ) of this equation will induce a local coordinate transformation for which 
H{x\ x 2 , x 3 ) has the form (A. 14), as desired. Note that F\ = H x i ^ 0, and hence F\ must have 
the same sign as e. ■ 

Local coordinates for which H has the form (A. 14) are determined up to transformations of 
the form (A. 5) with $> X 3 X 3 = 0, i.e., 

$(x 2 ,x 3 ) = $i(x 2 )x 3 + <f> Q (x 2 ). 
With <3> as above, the local coordinate transformation (A. 5) reduces to 

x 1 = x 1 + <Z> {x 2 ), x 2 = $i(x 2 ), x 3 = <S>' {x 2 ) + $' 1 (x 2 )x 3 . (A.16) 
With the assumption that H has the form 

#(x\x 2 ,x 3 ) = Fi(x 2 ,x 3 )x 1 + F (x 2 ,x 3 ), 
differentiating equation (A. 9) with respect to x 1 yields 

(*i)«» = Q 
V^Fi 

Therefore, 

Fi(x 2 ,x 3 ) =F 1 (x 2 ). 
Now equation (A. 8) reduces to 

F 1 (x 2 )x 1 + F (x 2 ,x 3 ) = (<I>;(x 2 )) 2 (F 1 (x 2 )x 1 + Fo(^^^ 3 ))+<I> '(x 2 )+$' 1 '(x 2 )x 3 , 

which, taking (A.16) into account, becomes 

F 1 (x 2 )x 1 + (Fx(x 2 )^ Q (x 2 ) + F (x 2 ,$ (x 2 ) + $' 1 (x 2 )x 3 )) 

= (^ 1 (x 2 )) 2 F 1 (^ 1 (x 2 ))x 1 + (^{x 2 )) 2 F (^ 2 )^ 3 ) + $oV) +&l{x 2 )x\ (A.17) 
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Equating the coefficients of x 1 on both sides yields 

F 1 (x 2 ) = (d>' 1 (x 2 )) 2 F 1 ($ 1 (x 2 )). 
Thus any solution $i(x 2 ) of the equation 



<S>[(x 2 ) = ^eF 1 (x2) 
will induce a local coordinate transformation for which 
Fi(x 2 ) = e = ±1. 

Local coordinates for which iq(x 2 ) = e are determined up to transformations of the form (A. 16) 
with 

for simplicity, we will assume that ^[(x 2 ) = 1. Then 

$ (x 2 , x 3 ) = x 2 x 3 + ax 3 + <J> (x 2 ) 
for some constant a. With <]? as above, the local coordinate transformation (A. 5) reduces to 

x 1 = x 1 + $o(x 2 ), x 2 = x 2 + a, x 3 = x 3 + & (x 2 ). (A. 18) 

Now equation (A. 9) takes the form 

(F ) x s = ec 3 . 
Therefore, 

F {x 2 ,x 3 ) =ec 3 x 3 + F 2 (x 2 ). 
Now equation (A. 17) reduces to 

e$o(x 2 ) + ec3$' (x 2 ) + F 2 {x 2 ) = F 2 (x 2 + a) + $o(x 2 ). 
Thus any solution $o( a;2 ) of the equation 

e<Z>o{x 2 ) + ec 3 $' (x 2 ) - $' '(x 2 ) = -F 2 (x 2 ) 
will induce a local coordinate transformation for which 

F 2 (x 2 ) = 0. 

Local coordinates for which F 2 (x 2 ) = are determined up to transformations of the form (A. 18) 
with 

e<D (x 2 ) +ec 3 ^ (x 2 ) -<J?' '(x 2 ) =0, 

i.e., 

$ (x 2 ) =b x e r ^ +b 2 e r * x \ 
where bi, b 2 are constants and 

n = £C3 + f^, r 2 = £C3 "f^. (A.19) 

Note that if e = 1, then n, r 2 are real and distinct; if e = —1, then n, r 2 may be real and 
distinct, real and equal, or a complex conjugate pair. 

To summarize, we have constructed local coordinates for which 

J(x 1 ,x 2 ,x 3 ) = ci, ^(x 1 ,^ 2 ,^ 3 ) = e(x 1 + c 3 x 3 ). 

These coordinates are determined up to transformations of the form 

x 1 = x 1 + bie rix2 +b 2 e r2X \ x 2 = x 2 + a, x 3 = x 3 + bme rix2 + b 2 r 2 e r2x2 . 
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A.2 c 2 = 

We will only give the details of the analysis for the case (A. 12); the case (A. 13) is similar. First 
we compute how the function (A. 12) transforms under a local coordinate transformation of the 
form (A. 5). 

Lemma A.2. There exists a local coordinate transformation of the form (A. 5) such that 
F Q (x 2 , x 3 ) = 0, i.e., 

x7(x\x 2 ,x 3 ) = Fi(x 2 ,x 3 )tan(c 3 5 1 ) + F 2 (x 2 ,x 3 ), (A.20) 
with Fx ^ 0. 

Proof. Substituting (A. 12) into the expression (A. 15) for .^(x 1 , x 2 , x 3 ) yields 
Fx tan (cgx 1 + F ) + F 2 



(F x tan(c 3 x 1 + F ) + F 2 ) - § X 2 X 2 



§ X 3 X 3 (Fx tan(c 3 x 1 + F ) + F 2 ) + (($^) 2 - Q^Q&ssj 

_ Fx sin(c 3 x 1 + F Q ) + (F 2 - $ x 2 x i) cos(c 3 x 1 + F ) 

$ x 3 x3 Fx sin(c 3 x 1 +F ) + {<& X 3 X 3 F 2 + (%2 S a) 2 - $ X 2 X 2 $ X 3 X 3 ) cos(c 3 x 1 +F y 

Now, define functions i?(x 2 ,x 3 ), 0(x 2 ,x 3 ) by the conditions that 

<& x 3 x 3Fi = i?sinO, & X 3 X 3F 2 + (& X 2 X 3) 2 — ^> :r 2 :r 2$ : j3 x 3 = i?cosO; 

in particular, we have 

= tan~ 



(A.21) 



Then the denominator of the right-hand side of (A.21) can be written as 

Rcos (c 3 x x + F - 0) = iicos (c^x 1 + $ - x 3 $ 53 ) + F - 0). 
Therefore, H(x , x 2 ,x 3 ) is a linear function of the quantity 

tan (c 3 (x 1 + $ - x 3 ^) + F - 0) , 
which implies that 

F = c 3 (<f> - x 3 <f> s3 ) + F - &. 

Keeping in mind that is a second-order differential operator in <3>, the condition Fq = 
is a second-order PDE for the function <3?(x 2 ,x 3 ). Any solution <l?(x 2 ,x 3 ) of this equation 
will induce a local coordinate transformation for which H( 3 ) has the form (A.20), as 

desired. ■ 

Local coordinates for which H has the form (A.20) are determined up to transformations of 
the form (A. 5) with <l>(x 2 ,x 3 ) a solution of the PDE 

c 3 ($ - x 3 ^) -0 = 0. (A.22) 

Unfortunately we cannot explicitly write down the general solution to this PDE; however, a sub- 
set of the solutions is given by the family 

$(x 2 ,x 3 ) = x 3 $ (z 2 ), 

where $o is an arbitrary function of x 2 with <3?q(x 2 ) ^ 0. 
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With the assumption that H has the form (A. 20), equation (A. 9) becomes (recalling that 
c 2 = 0) 

Fi((Fx) x3 - 2c 3 x 3 ) tan (^x 1 ) + 2F 1 (F 2 ) x s - F 2 {F 1 ) x z - (Fi) x 2 = 0. 
Therefore, since F\ ^ 0, 

(Fi) x3 - 2c 3 x 3 = 0, 2F 1 (F 2 ) x3 - F 2 (F 1 ) x3 - {F x ) x 2 = 0. (A.23) 
The first equation implies that 

F^x^x 3 ) =c 3 (x 3 ) 2 + F l0 (x 2 ) 

for some function i ? io(x 2 ). 

Computations similar to those in the previous case show that, under a local coordinate 
transformation (A. 5) with $ = x 3 $o( a;2 )) we have 

%(x 2 ) 

Thus any solution &q(x 2 ) of the equation 



will induce a local coordinate transformation for which either _Fio(x 2 ) = or Fio(i 2 ) = ±03. 
Denote this constant by c 4 , so that we now have 

Fi(5 2 ,5 3 ) = c 3 (x 3 ) 2 + c 4 . 
Finally, the second equation in (A.23) becomes 
2 (c 3 (x 3 ) 2 + a)(F 2 ) x3 - 2c 3 x 3 F 2 = 0, 
which implies that 

F 2 (x 2 ,x 3 ) = F 20 (x 2 )v / C3 (x3)2 + C4 

for some function -F2o(x 2 ). We conjecture that the remaining solutions of (A. 22) can be used to 
normalize the function F 2 q (x 2 ) , but unfortunately we have been unable to complete this step in 
the analysis. 

To summarize, we have constructed local coordinates for which 
j( 1 2 _ ci cos(c 3 a; 1 ) 

U I X . X . X I — , ■ 

Vc 3 (c 2 (a; 3 ) 2 + c 4 ) 
H(x 1 ,x 2 ,x 3 ) = (c 3 (x 3 ) 2 + c 4 ) tan [c 3 x x ) + F 20 (x 2 ) yjc 3 (x 3 ) 2 + c 4 . 
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